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The  new  Point  Analysis  (PA)  process  was  developed  to  correct  specific 
deficiencies  in  the  previous  process  (AFGWC  Point  Analysis  Documentation, 

1 Aug  1974) . 


These  deficiencies,  identified  by  the  3rd  Weather  Wing  Point  Analysis 
Conference  in  March  of  1976,  fall  into  three  general  categories.  First, 
the  moisture  modeling  technique  had  a deficiency  of  moisture  through  the 
tropopause.  Second,  there  was  a lack  of  observed  density  data  for  the 
80  km  to  250  km  region.  Finally,  the  production  process  equired  an  excessive 
number  of  manual  tasks. 


The  new  PA  process  employs  a new  tropopause  moisture  modeling  technique; 
includes,  when  available,  rocket  sonde  (ROCOB)  observations;  and  reduces  the 
manual  tasks  to  directing  and  quality  controlling  the  automated  process. 


The  work  described  here  is  primarily  the  result  of  the  efforts  of 
Capt  Paul  T.  Nipko,  Special  Projects  Branch  of  the  Production  Division, 
AFGWC.  It  is  based  on  the  work  of  his  predecessor,  Capt  Robert  D.  Abbey; 
and  includes  the  tropopause  moisture  modeling  technique  developed  by  Capt 
Lawrence  D.  Mendenhall,  USAFETAC. 
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1.  INTRODUCTION 


The  Air  Force  Global  Weather  Central  (AFGWC)  Point  Analysis  (PA) 
is  a process  which  produces  a description  of  the  atmospheric  conditions 
above  and  around  a given  point  on  the  earth.  The  final  product  includes 
a vertical  profile,  a gridded  cloud  depiction,  and  a site  weather  pseudo 
surface  observation . 

The  point  analyses  are  produced  by  both  AFGWC  and  USAF  Environmental 
Technical  Applications  Center  (USAFKTAC) . The  algorithms  used  by  both 
centralized  production  systems  are  identical,  however  the  completeness  of 
their  respective  data  bases  can  differ  significantly.  In  each  case  the 
final  product  represents  a blend  of  computer-processed  information  and  the 
judgement  of  an  experienced  weather  analyst. 

The  PA  process  is  triggered  by  a user  request  specifying  a time  and 
location.  Normally,  USAFETAC  will  process  requests  received  later  than 
18  hours  after  the  event  time.  AFGWC  is  charged  with  responding  immediately 
to  the  more  timely  requests.  This  limitation  is  driven  by  the  lifetime 
of  the  data  base  at  AFGWC.  Only  those  data  which  are  available  in  the 
data  base  at  the  time  of  request  are  included  in  the  AFGWC  analysis. 

AFGWC  response  time  for  a PA  request  within  an  18-hour  window  is  rarely 
more  than  four  hours  and  is  usually  much  less. 

2.  INPUT 


The  user  specifies  the  location  (latitude  and  longitude)  and  the  time 
for  which  a point  analysis  is  required. 

The  PA  process  uses  observed  data,  analyzed  data,  and  climatological  data. 
Observed  data  include  RAOB,  PIBAL  and  ROCOB  observations;  surface  observations; 
and  certain  parameters  which  define  solar  flux  and  geomagnetic  field  intensity. 
Analyzed  data  include  the  upper  air  analysis  fields  and  the  three  dimensional 
nephanaLysis  (3DNeph) . Climatological  data  from  the  U.S.  Standard  Atmospheric 
Supplements,  1066,  are  used  to  extend  the  vertical  profile  above  the  level 
of  observed  data. 

3.  CONSTRUCTION  OF  VERTICAL  PROFILE 

The  construction  of  the  vertical  profile  is  an  eleven  step  process. 

These  steps  are: 

(1)  identify  a suitable  raob  to  use  as  a skeletal  structure  for 
the  profile, 

(2)  obtain  values  from  the  analysis  fields  for  data  missing  at 
mandatory  levels  of  the  profile, 

(3)  complete  the  basic  profile  through  various  analysis  techniques 
(interpolation,  extrapolation,  regression) , 

(4)  add  PIBAL  and  ROCOB  data  to  the  basic  profile, 
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(5)  add  climatological  data  for  the  5 millibar  (mb)  level  to  the 
top  of  the  profile  (if  necessary  to  extend  it  to  100K  feet) , 

(6)  interpolate  meteorological  values  to  levels  at  500  foot  intervals, 

(7)  compute  density  for  each  500  foot  level, 

(8)  compute  absolute  humidity  for  each  500  foot  layer, 

(9)  compute  moisture  scale  height  for  each  500  foot  layer  below 
50K  feet, 

(10)  compute  precipitable  water  content  for  each  500  foot  layer  below 
100K  feet,  and 

(11)  extend  the  profile  above  100K  feet. 


3 . 1 Establishment  of  Skeletal  Structure  - RAOB  Data 

The  first  step  in  constructing  a vertical  profile  is  to  search  the 
data  base  for  a suitable  sounding.  A suitable  sounding  is  defined  as  one 
that  is:  (1)  within  15  hours  of  the  requested  time,  (2)  within  twice  a 
specified  radius  (usually  150nm) , (3)  characterized  by  at  least  four  values 

of  dewpoint  depression,  (4)  at  least  as  high  as  a specified  pressure  level 
(usually  400mb) , and  (5)  validated  by  data  base  management  procedures.  An 
example  of  the  data  available  in  a suitable;  RAOB  is  shown  in  Figure  1. 

The  decision  logic  used  to  select  the  best  sounding  among  those  which  are 
suitable  is  shown  in  Figure  2. 

If  no  suitable  sounding  is  identified,  the  data  in  the  upper  air 
analysis  fields  are  used  for  the  skeletal  structure,  unless  the  user  has 
specified  that  analysis  fields  should  not  be  used. 

3.2  Addition  of  Analysis  Fields  - Basic  Profile 


Analysis  fields  are  used  to  fill  missing  data  at  the  mandatory  levels 
of  the  skeletal  profile  (unless  othersise  specified) . Analysis 
data  are  available  at  grid  points  in  the  AFGWC  Northern  Hemispheric,  Southern 
Hemispheric  and  Tropical  grids.  Figures  3,  4 and  5 show  the  areas  of  the 
world  which  are  included  by  each  of  the  above  grids.  Figure  6 depicts  the 
pressure  levels  and  data  which  are  normally  available  in  the  analysis 
fields.  Note  that  the  70,  50,  30,  20  and  10  millibar  (mb)  levels  are  only 
available  for  the  Northern  Hemispheric  and  Tropical  grids,  and  only  for 
twelve-hourly  data  times  (00Z  and  12Z) . If  data  available  in  the  analysis 
fields  are  missing  in  the  skeletal  profile,  the  values  are  obtained  through 
bilinear  interpolation  of  the  analysis  fields.  The  combination  of  analysis 
fields  with  the  skeletal  profile  is  henceforth  referred  to  as  the  "basic 
profile” . 


3 . 3  Completion  of  the  Basic  Profile 

This  section  describes  the  mathematical  techniques  and  meteorological 
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methods  used  to  obtain  values  for  missing  data.  The  basic  profile  may 
have  missing  data  at  significant  RAOB  levels.  If  analysis  fields  are  not 
used,  data  may  also  be  missing  at  mandatory  RAOB  levels. 

The  surface  elevation  is  a critical  value  to  the  interpolation/extra- 
polation schemes  used  to  obtain  values  for  these  missing  data.  If,  for 
any  reason,  the  surface  elevation  is  missing,  and  the  value  can  not  be 
interpolated  from  the  AFGWC  Fixed  Terrain  Fields;  then  this  RAOB  is  reclas- 
sified as  unsuitable  and  another  search  is  conducted  for  a suitable  RAOB. 

Completion  of  the  basic  profile  is  a five  step  process: 

(1)  Determine  the  surface  pressure  if  missing. 

(2)  Determine  missing  temperature  values  throughout. 

(3)  Determine  missing  moisture  values  throughout. 

(4)  Determine  missing  height  values  throughout. 

(5)  Obtain  horizontal  wind  components. 

3.3.1  Surface  Pressure 

The  upper  air  levels  of  the  basic  profile  are  quality  controlled  to 
insure  that  pressure  is  filled.  If  the  surface  pressure  is  missing  and  an 
analysis  fields  level  exists  below  the  surface,  the  pressure  is  interpolated 
from  the  subterranean  level  to  the  surface.  Otherwise,  the  pressure  is 
extrapolated  from  the  lowest  upper  air  level  of  the  RAOB.  The  extrapolation 
formula  is  essentially  the  same  equation  as  that  used  to  obtain  height 
values  for  the  U.S.  Standard  Atmospheric  Supplements,  1966,  (appendix  A. 2) . 

3.3.2  Temperature  Throughout 

There  are  two  methods  to  obtain  a value  for  a missing  surface  temperature. 
If  a subterranean  analysis  level  exists,  the  surface  temperature  is  inter- 
polated between  the  subterranean  level  and  the  lowest  upper  level  with 
temperature  data.  The  interpolation  scheme  assumes  that  temperature  varies 
directly  with  the  natural  log  of  pressure.  This  is  analagous  to  drawing  a 
straight  line  between  the  temperature  values  at  two  pressure  levels  on  a 
thermodynamic  diagram.  If  a subterranean  analysis  level  is  not  available; 
then  the  product  is  assumed  to  be  constant  from  the  surface  to  the  lowest 
upper  level  is  solved  for  SRM,  then  the  equation  of  state  for  the  surface  is 
solved  for  temperature. 

If  temperature  data  do  not  extend  to  the  top  of  the  basic  profile,  the 
highest  temperature  data  are  extrapolated  to  the  top.  The  extrapolation 
scheme  uses  a lapse  rate  based  on  the  highest  temperature  data  and  the 
climatological  temperature  value  for  the  120,000  foot  level  (appendix  A. 3). 

Once  the  temperature  values  for  the  surface  and  the  profile  top  are 
known,  all  missing  temperatures  are  interpolated. 
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3.3.3  Moisture  and  Height  Throughout 

There  are  seven  steps  required  to  complete  the  moisture  profile: 


dewpoint. 


(1)  Calculate  mixing  ratio  from  pressure,  temperature  and 


(2)  Complete  mixing  ratio  to  highest  moisture  data. 

(3)  Calculate  missing  heights  throughout. 

(4)  Determine  the  tropopause  moisture. 

(5)  Set  a constant  mixing  ratio  at  and  above  the  hygropause. 

(6)  Interpolate  missing  mixing  ratios. 

(7)  Compute  dewpoints  throughout. 

3. 3. 3. 1 Calculate  mixing  ratio  from  P,  T,  T^ 

Moisture  data  in  a RAOB  and  in  the  analysis  fields  are  given  as  dewpoint 
depression.  Dewpoint  is  calculated  from  temperature  and  dewpoint  depression. 
Vapor  pressure  is  calculated  for  these  levels  from  pressure,  temperature  and 
dewpoint  (appendix  A. 4. 2).  Finally,  mixing  ratios  for  these  levels  can  be 
calculated  from  pressure  and  vapor  pressure  (appendix  A. 4. 3).  Moisture  data 
for  levels  with  a temperature  less  than  -40Oc  are  discarded  as  unreliable. 

3 . 3 . 3 . 2 Mixing  Ratio  to  Highest  Mixing  Ratio  Data 

If  the  surface  mixing  ratio  is  missing,  an  attempt  is  made  to  interpolate 
a value  using  the  lowest  upper  level  with  moisture  data  and  a subterranean 
analysis  level.  If  no  subterranean  analysis  level  exists,  the  mixing  ratio 
is  assumed  to  be  constant  from  the  surface  to  the  lowest  upper  level  with 
moisture  data.  Other  missing  mixing  ratios  (below  the  highest  level  with 
moisture  data)  are  obtained  by  interpolation  between  the  closest  levels 
with  moisture  data.  The  interpolation  formula  is  derived  from  the  assumption 
that  the  change  in  the  natural  log  of  mixing  ratio  is  inversely  proportional 
to  the  change  in  height  (appendix  A. 4). 

3. 3. 3. 3 Height  Throughout 

Missing  height  values  are  extrapolated  from  the  nearest  level  with  pressure, 
height  and  moisture  data.  The  extrapolation  formula  is  the  same  as  that 
used  to  compute  heights  for  the  U.S.  Standard  Atmospheric  Supplement,  1966 
(appendix  A. 2) . 

3. 3. 3. 4 Tropopause  Moisture 


The  tropopause  of  the  basic  profile  is  determined  in  accordance  with 
AWSM  105-124.  If  the  tropopause  is  undefined  the  new  moisture  modeling 
technique  cannot  be  used.  The  RAOB  must  be  reclassified  as  unsuitable  and 
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another  search  is  conducted  for  a suitable  RAOB.  If  the  tropopause  is 
defined,  the  next  step  is  to  interpolate  pressures  and  temperatures  for 
the  three  levels:  (1)  1KM  above  the  tropopause,  (2)  1 KM  below  the 

tropopause  , and  (3)  3 KM  below  the  tropopause.  Then,  calculate  the  lapse 

rates  between  the  tropopause  and  each  of  these  three  levels.  These  lapse 
rates,  the  tropopause  temperature  and  the  natural  log  of  the  tropopause 
pressure  are  the  predictors  used  by  a regression  equation  to  determine  the 
tropopause  dewpoint  (appendix  A. 4.1).  The  tropopause  mixing  ratio  is 
calculated  from  pressure,  temperature  and  dewpoint. 

3 . 3 . 3 . 5 Mixing  Ratio  Above  Tropopause 

The  hygropause  is  defined  to  be  the  level  above  which  mixing  ratio 
remains  constant.  The  dependent  data  used  to  derive  the  regression  eguation 
of  4. 3. 4. 3 indicate  that  the  hygropause  is  consistently  located  1KM  above  the 
tropopause  (USAFETAC  Report  7584) . A constant  mixing  ratio  of  3 parts  per 
million  (ppm)  is  set  for  levels  at  or  above  the  hygropause. 

3. 3. 3. 6 Mixing  ratio,  highest  data  to  tropopause 

The  mixing  ratios  are  interpolated  for  levels  between  the  highest 
moisture  data  and  the  tropopause  (appendix  A. 4) . 

3 . 3 . 3 . 7 Dewpoint  Throughout 


All  missing  dewpoints  are  calculated  from  pressure,  temperature  and 
mixing  ratio  (appendix  A.  5). 

3 . 4 Addition  of  PIBAL/ROCOB  Data 


PIBAL  and/or  ROCOB  data  associated  with  the  RAOB  are  added  to  the  basic 
profile.  Figure  7 is  an  example  of  PIBAL  data  typically  associated  with  a 
RAOB.  Figure  8 is  an  example  of  ROCOB  data.  The  addition  of  PIBAL/ROCOB 
data  is  accomplished  in  6 steps: 

(1)  Merge  the  PIBAL/ROCOB  levels  into  the  basic  profile. 

(2)  Combine  duplicated  levels. 

(3)  Interpolate  pressure  for  levels  added  below  the  basic  profile 

top. 

(4)  Interpolate  temperature  for  added  PIBAL  levels. 

(5)  Determine  moisture  values  at  added  levels. 

(6)  Calculate  pressure  for  ROCOB  levels  added  above  the  basic 
profile  top.  The  profile  resulting  from  these  six  steps  will  be  referred  to 
as  the  "merged  profile" . 
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3.4.1  Merge  the  P IBAL/ROCOB  Levels  Into  the  Basic  Profile 

PIBAL/ROCOB  levels  are  fitted  into  the  basic  profile  by  height.  Any 
pressure  data  which  are  included  in  the  PIBAL  or  ROCOB  are  discarded. 

This  is  necessary  to  assure  that  merging  the  different  data  types  does  not 
result  in  a pressure  inversion. 

3.4.2  Combine  Duplicated  Levels 

The  following  hierarchy  determines  which  values  are  used  if  two  or 
more  different  observation  types  report  the  same  physical  properties  at 
the  same  level: 

(1)  PIBAL  winds  are  used  over  both  RAOB  and  ROCOB  winds,  (2)  RAOB 
winds  are  used  over  ROCOB  winds,  and  (3)  RAOB  temperature  is  used 
over  ROCOB  temperature. 

3.4.3  Determine  pressure  for  Levels  Added  Below  the  Sasic 
profile  top 

Pressure  is  extrapolated  from  the  nearest  basic  profile  level  to  the 
added  level  (appendix  A. 2) 

3.4.4  Interpolated  Temperatures  for  Added  PIBAL  Levels 

Missing  temperatures  are  interpolated  as  in  section  3.3.2 

3.4.5  Determine  Moisture  Values  at  Added  Levels 


For  levels  added  below  the  basic  profile  top  missing  mixing  ratios  are 
interpolated  (appendix  A. 4) . ROCOB  levels  added  above  the  top  of  the 
basic  profile  are,  by  definition,  above  the  hvqropause.  Therefore,  a constant 
mixing  ratio  of  3ppm  is  used.  Missing  dewpoints  below  100,000  feet  are 
calculated  from  pressure,  temperature  and  mixing  ratio  (appendix  A. 5). 

No  attempt  is  made  to  extend  the  dewpoint  profile  above  100,000  feet. 

3.4.6  Calculate  Pressure  for  ROCOB  Levels  Added  Above  the  Basic 
profile  top 

The  gas  constant  for  moist  air  is  a function  of  mixing  ratio  alone. 
Therefore,  pressures  for  ROCOB  levels  (added  above  the  basic  profile  top) 
can  be  calculated  using  the  equation  of  state.  A check  is  required  to 
determine  if  the  pressure  at  the  top  of  the  basic  profile  is  less  than 
the  first  ROCOB  level  pressure  calculated  from  the  equation  of  state. 

If  a physical  inconsistency  exists  a new  pressure  value  is  extrapolated  for 
the  top  of  the  basic  profile. 

3.4.7  Wind  Components 

For  each  level  with  wind  data,  the  wind  direction  and  speed  are  resolved 
into  horizontal  components.  Values  for  each  wind  component  are  interpolated 
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to  those  levels  missing  wind  data.  A linear  interpolation  in  height  is 
used.  Wind  components  are  extrapolated  to  onLy  two  levels  above  the  highest 
data.  A linear  extrapolation  is  made  for  each  component,  using  the  average 
vertical  rate  of  change  for  the  u and  v components  over  the  highest  three 
levels. 

3. 5 Extension  to  100K  Feet 

The  merged  profile  must  extend  to  a height  of  at  least  100,000  feet 
above  ground  level  (AGL) . If  it  does  not  already  do  so,  the  5mb  level 
(always  above  100,000  feet  AGL)  is  added  to  the  top  of  the  merged  profile. 
Height  and  temperature  values  are  extrapolated  to  this  level  (appendices 
A. 6 and  A. 3).  The  mixing  ratio  is  set  at  3ppm. 

3. 6 Transformation  to  500  Foot  Levels 

The  merged  profile  is  transformed  from  non-linear  pressure  levels  to 
linear  height  levels.  Values  of  pressure,  horizontal  wind  components, 
temperature  and  mixing  ratio  are  interpolated  to  500  foot  levels.  A dewpoint 
value  is  calculated  for  each  500  foot  level  from  pressure,  temperature  and 
mixing  ratio  (appendix  A. 5). 

3 . 7 Computation  of  Density  for  Each  500  Foot  Level 

Two  methods  are  used  to  calculate  missing  density  values.  Below  the 
lowest  ROCOB  data,  density  values  are  calculated  from  the  equation  of  state. 
Above  the  lowest  ROCOB  data,  missing  density  values  are  linearly  interpolated 
in  height  from  ROCOB  data. 

3 . 8 Computation  of  Absolute  Humidity  for  Each  500  Soot  Layer . 

The  absolute  humidity  (water  vapor  density)  is  calculated  from  mixing 
ratio  and  density  (appendix  A. 7).  The  values  of  mixing  ratio  and  density  at 
each  500  foot  level  are  assumed  to  be  constant  for  the  layer  above. 

3 . 9 Moisture  Scale  Height 

The  moisture  scale  height,  msh,  is  defined: 

msh  = m 1_. 

V (1.1  - rh) 


where  m = the  mass  of  mater  vapor  and  rh  = relative  humidity  (decimal  value) . 

The  expression  J 1 )_  weights  the  amount  of  moisture  in  the  layer  such 

( 1.1  - rh  ) 

that  the  value  for  moisture  scale  height  increases  the  closer  the  moisture  is 
to  saturation.  1.1  is  used  in  the  denominator  instead  of  1.  so  that  moisture 
scale  height  will  not  be  undefined  for  the  saturated  case.  The  five  steps 
necessary  to  calculate  moisture  scale  height  are:  (1)  calculate  the  mass  of 
vapor  in  a column  of  air  having  a horizontal  cross  section  of  one  square 
meter  (mv  = p (thickness)  , (2)  calculate  the  saturation  vapor  pressure, 

es  from  the  temperature  (goff-gratch  equation  - appendix  A. 4. 2),  (3)  calculate 
the  saturating  ratio,  ws,  from  pressure  and  saturation  vapor  pressure 
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(appendix  A. 4 . 3) , 


calculate  the  relative  humidity  from  mixing  ratio 


and  saturation  mixing  ratio  (rh  = w/wg) , and  (5)  calculate  the  moisture 
scale  height  for  the  layer. 

3. 10  Computation  of  Precipitable  Water  Content,  (PWC) , for  Bach 
500  Foot  Level 

There  are  three  steps  to  compute  the  PWC  for  each  500  foot  layer.  First, 
overlay  the  transformed  profile  with  the  merged  profile.  Second,  calculate 
the  PWC  for  each  layer  of  the  overlayed  profiles  from  the  pressures  at  the 
boundary  levels,  the  mixing  ratios  at  the  boundary  levels,  and  the  gravitational 
acceleration  at  the  mean  height  of  the  layer  (appendix  A. 8) . Third,  add  the 
PWC  of  the  appropriate  layers  of  the  overlayed  profiles  to  get  the  PWC  of  the 
500  foot  layers  of  the  transformed  profile. 

3. 11  Extension  of  Profile  Above  100, 000  Feet 

For  points  in  the  Northern  Hemisphere,  the  pressure  values  are  extended 
to  200,000  feet.  Temperature  and  density  values  are  extended  to  400,000 
feet.  Values  of  these  parameters  are  determined  for  each  20,000  foot  level. 

For  points  in  the  Southern  Hemisphere,  no  extension  above  100,000  feet  is 
attempted. 

3.11.1  Extension  of  Pressure  to  200, 000  Feet 

If  ROCOB  data  art;  not  includeo,  pressures  are  extended  to  200,000  feet 
using  a hydrostatic  extrapolation  coupled  with  climatology  (appendix  B) . 

If  ROCOB  data  are  included  the  pressure  profile  will  normally  already  extend 
above  200,000  feet. 

3.11.2  Extension  of  Temperature  and  Density  to  400 , 000  Feet 

Temperature  and  density  values  (above  the  highest  data)  are  calculated 
using  climatological  parameters  derived  from  the  U.S.  Standard  Atmospheric 
Supplements,  1966  (appendix  C) . The  temperature  and  density  profiles  join 
smoothly  with  the  Jacchia  Model  boundary  conditions  at  120  KM  (Part  3, 
Atmospheric  Models  Above  120  KM,  U.S.  Standard  Atmospheric  Supplements,  1966). 

3.12  Summary  of  the  Vertical  Profile 

Only  specific  levels  from  the  completed  profile  are  included  in  the 
final  Point  Analysis  output  (Figure  9b) . For  points  in  the  Northern  Hemisphere, 
the  following  physical  * -operties  are  included:  (1)  pressure  up  to  200,000 
feet  AGL.  (2)  horizontal  wind  components  (highest  level  dependent  upon 
available  data),  (3)  temperature  up  to  400,000  feet  AGL,  (4)  absolute 
humidity  up  to  100,000  feet  AGL,  (5)  density  up  to  400,000  feet  AGL,  and 
(6)  precipitable  water  content  for  each  layer  below  100,000  feet  AGL.  For 
points  in  the  Southern  Hemisphere  the  same  physical  properties  are  included 
up  to  100,000  feet  AGL. 
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4.  CLOUD  DEPICTION 


A Cloud  Depiction  (CD)  is  provided  for  the  area  surrounding  the  input 
point.  The  CD  displays  values  of  low,  middle,  high  and  total  cloud  amounts 
in  eighths  of  coverage.  These  values  are  displayed  at  the  grid  points  of 
a polar  stereographic  projection  (data  grid) , displayed  for  a map  scale 
of  1:1,  875,000.  The  data  grid  is  true  at  60°N  or  60°S  latitude  (where 
the  data  grid  spacing  is  approximately  25  ran) . The  display  is  centered  on 
the  input  point  and  covers  a square  area  with  sides  equal  to  six  times  the 
grid  spacing  at  the  input  point.  A reference  grid  (equally  spaced  in  earth 
distance)  is  superimposed  on  the  data  grid.  The  reference  grid  is  always 
100  nautical  miles  on  a side  with  grid  spacing  of  25  nautical  miles.  It  is 
centered  on  the  input  point  and  squared  to  the  printed  page.  In  the  northern 
hemisphere,  north  is  towards  the  top  of  the  grid.  In  the  southern  hemisphere , 
south  is  at^  the  top  of  the  grid.  The  cloud  coverage  amounts  are  plotted  top  to 
bottom;  high  through  low,  respectively.  The  total  cloud  amount  is  plotted  to  the 
right  of  the  middle  cloud  amount  (figure  9a) . 

4 . 1 Source  Data  for  the  CD 

Data  for  the  CD  are  extracted  from  the  Three  Dimensional  Nephanalysis 
( 3DNEPH)  portion  of  the  AFGWC  data  base.  The  3DNEPH  (AFGWC  TM  70-9  and 
AFGWC  TM  71-2)  is  generated  every  three  hours  (beginning  with  00Z) . It  is 
produced  from  a combination  of  meteorological  satellite  observations  and 
conventionally  sensed  data.  The  3DNEPH  provides  percent  of  cloud  cover  for 
15  vertical  layers  of  variable  thickness.  The  lower  six  layers  are  terrain 
following,  while  the  upper  nine  layers  are  referenced  to  mean  sea  level  (MSL) . 

The  actual  layers  are: 

(1)  0 - 150  FT  AGL, 

(2)  150  - 300  FT  AGL, 

(3)  300  - 600  FT  AGL, 

(4)  600  - 1000  (1M)  FT  AGL, 

(5)  1M  - 2M  FT  AGL, 

(6)  2M  - 3.5M  FT  AGL, 

(7)  3. 5M  - 5M  FT  MSL, 

(8)  5M  - 6. 5M  FT  MSL, 

(9)  6.5M  - 10M  FT  MSL, 

(10)  10M  - 14M  FT  MSL, 

(11)  14M  - 18M  FT  MSL, 

(12)  18M  - 22M  FT  MSL, 
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(13)  22M  - 26M  FT  MSL, 

(14)  26M  - 35M  FT  MSL,  and 

(15)  35M  FT  MSL  and  above. 

The  3DNEPH  also  provides  a total  percent  of  cloud  coverage  for  all  layers. 

4 . 2 Combination  of  3DNEPH  Layers  to  Determine  Low,  Middle  and  High 
Gloud  amounts 

The  15  layers  of  the  3DNEPH  are  combined  to  determine  the  amounts  of 
low,  middle  and  high  clouds.  The  limits  of  these  three  cloud  categories 
vary  somewhat  with  surface  elevation.  The  approximate  limits  are:  (1)  0 to 
6000  feet  AGL  for  low  clouds,  (2)  6000  to  20000  feet  AGL  for  middle  clouds, 
and  (3)  above  20,000  feet  AGL  for  high  clouds.  The  combining  process  is  an 
approximation  which  makes  use  of  a "stacking"  factor  (appendix  A. 10).  This 
factor  is  a measure  of  the  degree  to  which  individual  cloud  layers  are  vertically 
separated. 

5.  SURFACE  OBSERVATIONS 

Data  from  the  four  closed:  surface  observations  may  be  tabulated  in  paragraph 
G.  of  the  final  product.  The  maximum  distance  from  the  site  to  any  tabulated 
surface  observation  is  approximately  375  nm.  Figure  9c  has  an  example  of 
such  a table. 

6.  SITE  WEATHER  PSEUDO  OBSERVATION 


A pseudo  surface  observation  is  created  by  using  data  available  from 
the  3DNEPH  and  the  four  closest  surface  observations  described  in  section  5. 
(See  also  8.2) 

6. 1 Surface  Visibility 

The  surface  visibility  is  a distance  weighted  average  of  the  visibilities 
reported  at  the  four  closest  stations.  If  no  surface  observation  data  is 
available,  the  visibility  is  reported  as  99.99.  Visibility  is  in  nautical 
miles. 


6 . 2 Cloud  Types  (low,  middle,  high) 

The  cloud  types  present  in  the  3DNEPH  data  base  at  the  closest  grid  point 
are  assumed  to  be  present  at  the  input  site. 
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6 • 3 Cloud  c°ver  (in  eighth^  for  Each  cloud 
(low,  middle,  high ) 


C a teg 


The  low,  middle  and  high  cloud  amounts  calculated  in  section  5 (for 
the  nearest  3DNEPH  grid  point)  are  assumed  to  be  valid  for  the  input  site. 


6 . 4  Layer  Limits  (minimum  bases  and  maximum  tops) 


The  3DNEPH  layers  for  the  nearest  3DNEPH  grid  point  are  examined  to 
determine  which  layers  make  a non  zero  contribution  to  the  cloud  amounts  for  each 
cloud  category  (low,  middle,  high) . The  extremes  of  the  highest  and  lowest 
layers  used  for  each  category  are  used  as  the  layer  limits  for  that  category, 
the  heights  are  in  hundreds  of  feet. 


6. 5 Present  Weather  (for  the  input  time) 

The  weather  reported  at  the  nearest  surface  reporting  station  is  assumed 
to  be  the  weather  at  the  site. 

6.6  Surface  Wind 

Wind  direction,  speed,  and  gust  speed  are  distance  weighted  averages 
of  the  four  closest  surface  stations.  If  no  wind  data  is  available,  direction 
and  speed  will  both  be  "999".  Missing  gust  data  or  no  gusts  are  both  coded 
as  "***".  wind  direction  is  in  degrees.  Wind  speed  and  gust  speed  are  in 
knots. 


7.  FINAL  PRODUCT 


The  final  product  contains  seven  paragraphs. 

7 . 1 Paragraphs  A through  F . 

Paragraphs  A through  F contain  the  following t (A)  -either  a six  character 
site  code  or  the  latitude  and  longitude;  (B)  -either  blank  or  the  input  time 
of  interest;  (C)  -cloud  depict; on  discussed  in  section  4;  (D)  -the  site  weather 
pseudo  observation;  (E)  the  precipatible  water  information  from  the  vertical 
profile  discussed  in  section  3;  (F)  -all  the  remaining  parameters  from  the 

vertical  profile. 

7. 2 Paragraph  G. 

Paragraph  G begins  with  a number  of  plain  language  automated  remarks. 

These  remarks  describe  the  types  and  times  of  the  data  used  to  create  the 
vertical  profile.  If  a RAOB  was  used,  the  distance  and  bearing  to  the  RAOB  will 
appear.  The  RAOB  station  number  may  also  appear.  The  last  automated  remark 
has  the  characters  "MARK"  in  columns  75  through  78.  This  line  may  be  followed 
by  the  surface  observations  table.  The  surface  observation  table  can  be 
further  controlled  to  display  or  not  display  station  identification  information. 
The  next  record  to  appear  is  the  solar  flux  and  geomagnetic  field  intensity 
parameters.  Then  comes  manual  forecaster  remarks  in  plain  launguage.  These 
manual  remarks  may  discuss  the  representativeness  of  the  RAOB  or  any  other 
information  the  forecaster  wants  the  PA  customer  to  know.  The  last  record  of 
the  PA  gives  the  vertical-profile-software  version  number. 
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8.0  FORECASTER  MODIFICATIONS 


There  are  two  basic  categories  of  action  the  forecaster  can  take  to 
control  the  production  process. 

8 . 1 Control  of  Automated  Process 

The  options  which  the  forecaster  can  exercise  to  control  the  automated 
process  are: 

(1)  Remove  a maximum  of  five  upper  air  stations  from  consideration 
as  a suitable  RAOB. 

(2)  Specify  the  upper  air  station  to  use  as  the  suitable  RAOB. 

(3)  Change  the  Size  of  the  radius  within  which  the  search  for  a suitable 
RAOB  is  conducted. 

(4)  Change  the  pressure  level  above  which  a suitable  RAOB  must  extend. 

(5)  Specify  that  analysis  fields  are  not  to  be  used. 

(6)  Specify  that  only  analysis  fields  are  to  be  used. 

8 . 2 Manual  Changes 

The  forecaster  can  adjust  the  site  weather  pseudo  surface  observation 
in  any  way  he  feels  necessary.  He  may  also  insert,  in  paragraph  G.,  any 
plain  language  remarks  that  he  feels  may  help  the  customer  in  using  the  PA 
product. 

9.  SUMMARY 

The  PA  is  an  approximation  of  the  atmospheric  structure  and  the  weather 
phenomena  surrounding  a given  point.  Exactly  how  well  the  PA  represents  truth 
is  unknown.  However,  it  is  the  best  estimate  available.  The  total  PA  process 
makes  optimum  use  of  data  automation;  while  providing  the  capability  for  a 
trained  weather  analyst  to  control  and/or  modify  the  automated  process.  The 
result  is  a final  product  of  higher  quality  than  either  the  analyst  or  the 
computer  could  produce  separately. 
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Figure  B.  This  is  an  example  of  useable  Rocob  data. 
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FIGURE  9A:  Paragraphs  A.  - D.  of  an  example  PA. 
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Figure  9C.  Paragraph  G of  an  example  PA. 
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The  theoretical  development  in  this  appendix  is  confined  to  formulas 
utilized  below  100,000  feet.  Formulas  used  to  obtain  values  above  100,000 
feet  are  discussed  in  Appendices  B and  C. 

A. 1 Definitions 


A = g auto- convective  temperature  gradient,  the  temperature 

gradient  of  a constant  density  atmosphere  (°K.cm”l) 


direction  form  which  wind  is  blowing  (degrees 
measured  from  north) 


e 


partial  pressure  due  to  water  vapor 


"v 

£=  Rd  = 0.62197  . . . 

l 


saturation  vapor  pressure  at  the  icepoint  temperature 
(p=6. 1071  mb) 

saturation  vapor  pressure  (100%  relative  humidity) 

saturation  vapor  pressure  at  the  steampoint  temperature 
(p=1031 . 246  mb) 

ratio  of  dry  air  gas  contant  to  water  vapor  gas 
constant 

latitude 


2 

g acceleration  due  to  gravity  (cm/sec  ) (a  function 

of  latitude  and  height) 

G = = 3.085462(10)_6+2.27(10)-9  cos  (2I)"2  (10)"12  cos  (41) 

gravity  gradient  at  z=0  (List  1963) 


g({)  = 978. 0356(1. +0.0052885  sin2  £-0. 0000059sin2 (2$)  

g at  latitude  Z=0 

(m.sec-^)  (List  1963) 

^ Tv  vertical  temperature  gradient  in  the  direction 

^ z of  a decreasing  temperature  (°K/km) 
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height  scaling  parameter  (cm) 


latent  heat  (erg.gm  ^ or  cal  gm  *) 


Molecular  weight  of  dry  air 
(28.96(10-^)  Kg/mole  below  80Km) 


. atmospheric  pressure  (mb) 


precipitable  water  content  (CM) 


Q = w_(100) relative  humidity  (percent ) 

ws 


ViL 

e 

c 


specific  humidity  (gm  H^O  vapor/Kg  moist  air) 
total  atmospheric  density  (gm.cm-^) 


dry  air  density  (gm.cm-^) 


water  vapor  density  or  absolute  humidity  (gm.cm-3) 


r = 2g ($) 
G 


effective  earth  radius 


gas  constant  for  dry  air  (2.8704 (10)^  erg.gm_^.°K 


gas  constant  for  moist  air 


gas  constant  for  H2O  vapor  (4. 6150  (10) ^ erg.gm-'*' .°K  ^) 


absolute  temperature  (°K=Ti+aC) 


dewpoint,  the  temperature  to  which  moist  air  must 
be  cooled  at  constant  pressure  to  reach  saturation 
(°K) 


ice-point  temperture  (273.16°K  or  0°C) 
dewpoint  depression  (°K) 
steam-point 


T = T-Tj 
p d 


H=  "-2»I  + + t ri 

(rr('r+*)  G-f(r**) 


Z = rH 


gQr  - h 


where  Z is  geometric  altitude. 
The  perfect  gas  law  yields: 


(where  M is  the  molecular  weight  of  dry  air  = 28.96(10-3)K£  below  80Km.) 

mole 

While  the  hydrostatic  equation  can  be  expressed: 


therefore, 


Ip  = - f gdz, 


d p = -_M P g dz . 


Once  again,  substituting  for  g yields. 


dp  = - mp  gQr2  dz 

R*t7  (r+z)2 


rv  dp  = - Mgc 


(r+Z) 2 
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Assuming  T for  a layer  is  a good  approximation  to  Tv,  we  integrate  over 
the  layer  to  get: 


R*T 


r't/  y dp:  ~ r^/( 

f(  1 *-  0 


° (V'+i,')  (r+i,) 

+ 012^  _ /Wj, 
(r+^  (r+i.) 


**tA»(S*  -&f^[ 


Substituting  (2)  into  the  equation  above  yields 


R*T  In 

Let  H2  - = h,  then 

In  ( f 2/  Px)  = 

or 


r / P j)  - - GpH  Wj-Hj) 

-G-r^VR«T 

( P 2/  P l)  = exp  (rp  ^ Vl  / R*  T ) 
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(5) 


(6) 


k 


(6)  can  be  rewritten  as 


If  we  approximate  G with  g 


* ' 


P»/P,  = e*rf  . 


which 


is  identically  equation  (1.13)  from  the  U.S.  STANDARD  ATMOSPHERIC  SUPPLEMENTS, 
1966.  Equation  (6)  is  used  to  extrapolate  missing  surface  pressure. 

Equation  (6)  may  be  written: 


G,M 


(7) 


(8) 


(7)  or  (8)  is  used  to  extrapolate  geopotential  height  to  a level  where  height 
is  missing.  Equation  (3)  is  used  to  convert  the  extrapolated  geopotential 
height  to  geometric  height. 

A. 3 Temperature  Extrapolation 

let:  climo  temp  @ 120,000  ft  = T^2q 

height  in  meters  of  120,000  ft  = Z 

height  in  last  temp  level  = ZT 

Li 

temperature  at  last  temp  level  = Tt 

Li 

pressure  at  last  temp  level  = PL 
pressure  at  120,000  ft  = P 

Equation  of  State  for  Moist  Air  is:  p _ eRATv  (9) 

Rewriting  (9)  and  substituting  into  (4)  for  £ gives: 


- - P 


1 


¥ - 
P RaTv 


(10) 
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Now  assuming  g and  constant,  integrate  (10)  to  obtain: 

Then  take  partial  of  (11)  with  respect  to  lnp  to  obtain: 

(12) 

Now  assume  iT/>(iv,|»)  = constant  over  interval  Z^  to  Z (this  is 

analagous  to  drawing  a straight  line  on  a thermodynamic  die^9am)j  and 


T-Tu  _ }T_ 


or  T = In  (P/P.)  ^ T 

<KM 

but  from  (12) : 


+ Tt 


(13) 


JI  . -IWTv  yr 
% & 


(14) 


However,  (11),  (12)  and  therefore,  (14)  are  based  on  the  assumption  that  g 
and  T are  constant  throughout  the  layer.  If  we  assume  g at  mean  layer 
height  is  representative  of  the  entire  layer,  and  assume  ?hat  (T  + T, „„)/2 


represents  for  the  entire  layer,  then 


120 


II  _ = -fcj  (Tu#+V) 


(15) 


Now  assume 


11 

& 


at 

A2- 
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Then 
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Xuf) 


U 

a8»* 


(16) 


Substituting  (16)  into  (13)  yields: 


Kj  (Tt*-T,?e) 


A. 4 MIXING  RATIO  INTERPOLATION: 


(17) 


Assume  the  change  in  natural  log  of  mixing  ratio  is  inversely  proportional 
to  change  in  height;  then 


CJ(i*.w)=^ 


(18) 


It  has  been  shown  that  (12) 

Ofy  ^ 

Assumming  ^/Tv  constant  yields, 

t[l(*-p)l= 

Combining  (18)  and  (19)  obtains 


c»  p) 


A lv\w  __  y^  (20) 

A ^ p 


Let  subscript,  t,  denote  the  upper  boundry  level  for  a layer  and  subscript, 
b,  denote  the  lower  boundary  level  of  the  same  layer. 


Then,  A lnw  lnw  - lnw  ln(w  /w.) 

z b t b 

^lnp  lnpfc  - lnpb.  In  (Pt/Pb) 
Now  for  any  level  in  the  layer 

>»(w/wt) 

*«(!*/ po  * 

or  w=  ( P/ptf 


(21) 


A. 4.1  Tropopause  Dewpoint 

o 

Let:  T = Tropopause  Dew  Point,  C. 

G1A  = Average  Lapse  Rate  for  lKm  above  the  tropopause. 

GIB  = Average  Lapse  Rate  for  lKm  below  the  tropopause. 

G3B  = Average  Lapse  Rate  for  3Km  below  the  tropopause. 

Tfc  = Tropopause  temperature,  °C. 

PLOG^  = Natural  Log  of  tropopause  pressure  in  mb. 

The  regression  equation  used  to  fill  the  tropopause  dewpoint  is 

Tt  = Rt  + CT  * G1A  + E * GIB  + F*G3B  + R * Tfc  + S * PLOGfc. 
where  the  constants  are: 

RI  = -88.05397 
T = 0.45038 

E = - 0.84287 

F = - 1.1  8564 

R = 0.41794 

S = 6.384336 

* = multiplication 

For  information  concerning  the  development  of  this  regression  equation  see 
USAF  ETAC  Report  7584,  20  Aug  1975  and  USAF  ETAC  Report  8035,  Aug  1976. 
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A. 4. 2 Calculate  Saturation  Vapor  Pressure 

Computations  based  on  Gof f-c.ratch  equations  for  saturation  vapor 
pressure  over  ice  or  water  - see  Smithsonian  METEOROLOGICAL  TABLES,  LIST, 
6th  Edition,  page  350. 

Over  water  the  vapor  pressure  equation  is: 

e = 10LOG10W  ( 373 . 16/Td)  (2; 


where  Td  is  in  °K,  e is  in  mb,  and  LOG10W  (F)  =-7.90298  (F-l)  + 5.02808  LOG10F 

-1 . 3816  ( 10-7 ) (3.0  (11  * 344  (1-1/F)  )_X) 

+8.  1328 (10~3)  (io(3'4914'J(1-F)  *-1) 

LOG  (1013.246). 

+ 10  (23) 

Over  ice  the  vapor  pressure  equation  is: 

es  = 10L°Glct273-  16//Td)  , (24) 

Where  LOG10I  (F)  = -9. 09718 (F-l)  -3.56654  LOG-^F) 

+. 876793 ( 1-1/F ) + L0G10 (6. 1071) . (25) 

Vapor  Pressure  is  calculated  over  water  for  temperatures  above  -40°  to  be 
consistent  with  raob  measurements. 

A. 4. 3 Calculate  Mixing  Ratio 

The  equation  of  state  of  water  vapor  is  e = ^vRvT.  (26) 

The  equation  of  state  of  dry  air  is  = ^!RdT  (27) 

The  total  pressure  is  the  sum  of  partial  pressures 

p = e + pd,  (28) 

or  p-e  = pd. 

Dividing  (26)  by  (27)  yields 

e/Pd  = fv^T 

?dRdT  ' 

or  e/pd  = <{v)  (Rv) 

"6  *d  • (29) 
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Let  M = Mass  of  vapor  for  volume  V, 
v 

= Mass  of  dry  air  for  volume  V. 

Then  O = M /V, 

% v v 

and  fa = Vv- 

Therefore,  (29)  becomes 
e/Pd  = (Mv/V)  ( 1 ) 


(Md/V)  Rd/Rv 


but  M /M  = W, 
v a 

and  P = P-e. 
d 

If  Rd/Rv  = GASRAT, 

then  (30)  becomes 

e (w)  ( 1 ) 

P-e  GASRAT  , 

or  (e)  (GASRAT)  =W 
(p-e) 

A. 5 Calculate  Dewpoint 


R = ( M 
m v 


) R + ( M, 


) 


M +Mj 
v d 


Mv  + Md 

(Haltiner  and  Martin  (1957),  Page  24.) 


where  R = gas  constant  for  moist  air, 
m 

R,  = gas  constant  for  dry  air, 


Rv  = gas  constant  for  water  vapor, 
Mv  = mass  of  water  vapor, 


Md  = mass  of  dry  air. 

Multiplying  each  term  of  (32)  by/ 1/M^yields 

\ 1/M , 


ir 


^ = (VMd>  *v  + 


) 


Mv/Md+1 


Mv/Md+1 


(30) 


(31) 


(32) 
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r 


Since  w = M /M, , 
v a 


Rm  = ( w ) Rv  + ( 1 ) R 


d ' 


w + 1 w + 1 

The  Equation  of  State  of  Moist  Air  is  p = ^sRmT' 
or  £-  RmT/ p • 

For  any  volume  (V)  w can  be  alternately  defined  as 


w T—\  - e* 

\Md /v/  ~ 7* 


(33) 

(34) 

(35) 

(36) 


Since  the  total  density  equals  the  sum  of  the  densities,  we  have: 


Therefore,  (36)  yields 


(37) 

then  (26)  yields  vapor  pressure,  e. 

Saturation  vapor  pressure  (eg)  is  calculated  from  (22)  or  (24)  using  the 
temperature  T as  a first  approximation  to  T , . The  equation  which  relates 
these  two  values  is  the  Clausius-Clapeyron  Equation  (LIST,  1963) 

1 de  L 
s 

dT*“R^T2.  (38)  (38) 

By  approximating  de  /dT  with  A e / At  and  using  a value  for  latent  heat, 
(L)  based  on  the  temperature,  (3§)  becomes 


AT = 


RvTx4e> 


(39) 
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which  can  be  used  to  extrapolate  the  missing  height  Z 2_ 

The  final  two  steps,  (43)  and  (44),  are  repeated  until  the  change  in 
height  is  less  than  10-5  meter. 

A. 6. 2 Method  2 

If  TTLAP  in  (44)  is  zero,  method  2 is  used  to  extrapolate  heights. 

It  has  been  previously  shown  that  (12) 
j(ln  p)  = -g  JZ 


Using  the  delta  approximation  yields 

^ Inp  = -g  &Z 

RdT 
a v 

Approximating  g for  the  layer  by  g at  p-^  and  T for  the  layer  by  T 
at  p-^  we  have: 

4 


<?>r 


Tv,  j^CFVr,) 


(45) 


The  value  of  Z obtained  from  (45)  is  used  to  calculate  a new  value  for  g 
and  equation  (43)  is  solved  for  TTLAP.  Then  method  1 is  used  to  obtain 
a new  value  of 

A. 7 Absolute  Humidity 


The  definition  of  mixing  ratio  is 


w= 


(46) 
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Total  density  equals  the  sum  of  the  densities  of  the  constituents. 

Equation  (47)  may  be  rewritten 


ej  = e-ev 


Substituting  (48)  into  (46)  yields 


e. 


(l+w) 

However,  is  specific  humidity,  . Therefore, 

<te=ev 

and,  averaged  over  a layer. 


?e=ev 


Then 
we  have 


letting  q = (eTteRy 


(47) 


(48) 
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A. 8 Computation  of  Preciptable  Water  Content 

Let  average  density  of  water  vapor  in  a column  of  air  through  a given 

layer  (gm/cm^). 

Assume  a horizontal  cross  section  of  ) cm2  for  the  column. 


rnv=  Z,  a?  (1) 


Where,  = mass  of  water  vapor  in  the  column  of  air  for  this  layer  (gm) , and 
z = the  thickness  of  the  layer  (cm) . 


?v-(Cvt  +ev6)/a 


If  the  water  vapor  in  the  column  (niy)  is  condensed  into  a vessel  of 
equal  cross  section,  we  know  that 

mv  = C W ^^1,0*1 

Where,  “ density  of  liquid  water  (lgm/cm^),  and 

height  of  the  liquid  water  in  the  vessel  (precipitable  water 

content) . 

Use  PWC  for  then  combining  (50)  and  (51)  yields; 

ev  Air-  fVVC  (52 


pwc  = 

Cw 


However,  integrating  (4)  yields 


Substituting  (54)  into  (53)  yields 


£y  f —A  P 
ew  ^ l $ ' 
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A . 9 Interpolation  for  Site  Pseudo  Surface  Observation  Visibility  and  Wind 

Interpolation  from  the  four  surrounding  surface  observations  is  done 
by  a distance  weighted  average.  The  weight  factor  is 

F*  «*F  [ ” Cat)*”], 


where  F is  the  weight  factor  and  d is  the  distance  from  the  station  to 
the  point.  Notice  that  if  a station  is  at  the  input  point,  d would  be 
zero  and  F would  be  one.  As  d becomes  larger  F decreases,  with  zero 
as  a limit.  Distance  is  divided  by  25  to  keep  F significantly  larger 
than  zero  for  the  expected  range  of  distance  values. 

Visibility  is  calculated  with  the  following  equation: 

^ I 

' l F i > 

i=l 

where  V is  the  calculated  visibility,  V.  is  the  visibility  reported  by  a 
station,  and  F.  is  the  weight  factor  for  a station.  Variable  n can  vary  from 
1 to  4,  depending  on  how  many  surface  reports  are  tabulated.  If  there  is  no  data 
available,  visibility  is  set  at  99.99. 

The  wind  reported  by  each  surface  station  is  broken  into  east-west  components 
(ui)  and  north-south  components  (v^) . Then  a weighted  average  is  determined 
for  both. 


These  average  components  are  used  to  recreate  a final  direction  and  speed. 

If  wind  gust  data  is  available,  average  gust  is  determined  using  an  analogous 
equation.  However,  gust  data  is  reported  and  available  only  for  North 
America. 


Users  should  be  aware  that  the  surface  wind  in  paragraphs  D and  F are 
determined  in  completely  different  ways.  Therefore,  significant  differences 
should  be  expected. 

If  no  wind  data  is  available,  the  direction  and  speed  both  become  999. 
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A . 1 0 Low,  Middle,  and  High  Cloud  Amounts 

There  are  15  layers  in  the  3DNEPH  cloud  analysis.  These  layers  are  grouped 
into  low,  middle,  and  high  layer  groups.  Then  the  layers  within  a group  are 
combined.  First,  the  two  lowest  layers  in  a group  are  combined.  Then,  the 
next  layer  up  is  combined  into  the  result  of  the  first  combination,  and  so 
on  until  all  the  layers  of  a group  are  combined  into  a final  result  for  the 
group. 

The  process  of  combining  two  layers,  or  of  combining  a combination  with  the 
next  layer  up,  uses  a dimensionless  stacking  factor,  R.  R equal  to  zero 
indicates  that  the  clouds  in  the  next  layer  up  are  all  stacked  directly 
above  the  clouds  below.  R equal  to  one  indicates  that  the  clouds  in  the  next 
layer  up  are  totally  random  in  relation  to  the  clouds  below. 

The  difference  in  height  between  two  layers  (D)  is  used  to  determine  R.  If 
D is  less  than  610  meters,  R is  set  at  zero.  For  D greater  than  610  meters, 

R = D-610 
6100 

The  equation  used  to  determine  the  total  cloud  coverage  for  the  combination 
of  any  two  layers  is 

T = MAX  + (1-MAX) (MIN) R. 

where  T is  the  total  coverage,  MAX  is  the  larger  of  the  two  coverages,  and 
MIN  is  the  smaller  of  the  two  coverages. 
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APPENDIX  B 


EXTENSION  OF  PRESSURE  TO  200,000  FEET 
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The  hydrostatic  equation  can  be  expressed  by 


where  C is  a constant  of  integration  and  the  other  symbols  are  defined 
in  Appendix  A.  Changing  the  logarithm  base  from  natural  to  common 
yields 


log 


10 


P = Az  + B, 


w^trc.  A = ~ lojj<0«  j Av*df  C )o^loe 


Given  values  for  A,  which  is 

essentially  the  inverse  of  the  scale  height,  p can  be  extrapolated  to 
any  height. 

The  technique  programmed  for  the  PA  process  computes  B for  the  pressure 
at  100,000  feet  and  then  computes  the  pressure  for  each  20,000  feet  to 
200,000  feet.  Values  for  A,  taken  from  the  COSPAR  International  Reference 
Atmosphere  (1965)  for  60°  North  latitude  are  tabulated  below. 


Month 


A (X10  3) 


Jan 

-1.924 

Feb 

-1.918 

Mar 

-1.874 

Apr 

-1.789 

May 

-1. 705 

Jun 

-1.651 

Jul 

-1.645 

Aug 

-1.690 

Sep 

-1.733 

Oct 

-1.800 

Nov 

-1.847 

Dec 

-1.896 

appendix  c 


EXTENSION  OF  DENSITY  AND 
TEMPERATURE  TO  400,000  FEET 
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This  appendix  presents  a method  for  extending  the  temperature  and 
density  profiles  from  100,000  feet  to  400,000  feet,  using  climatological 
parameters. 

To  extend  the  density  profiles  from  100,000  feet  to  400,000  feet,  mean 
scale  heights  were  derived  from  the  U.S.  Standard  Atmospheric  Supple- 
ments (1966),  the  equation  used  is: 

(1) 

f (Z)  = f(Zo)exp^-  -tz.  ~ £&> 
where  ’’ 

f(zo)  is  the  density  at  height  zQ, 

^ (Z)  is  the  density  at  height  Z (Z>ZQ),  and 
H is  the  scale  height. 

The  density  profiles  are  broken  into  layers  for  which  (1)  is  valid 
(FiguresCl  andC2) . Table  C 1 presents  the  calculated  scale  heights  for 
the  appropriate  layers  for  given  months  and  latitudes.  The  error  in 
reproducing  the  densities  from  30-120  Km  in  the  U.S.  Standard  Atmospheric 
Supplements  (1966)  by  using  (1)  and  the  scale  heights  in  Tabled,  is  less 
than  10  percent. 

There  is  contradictory  evidence  concerning  the  existence  of  an  isopycnic 
(constant  density)  layer  in  the  vicinity  of  90  Km.  For  this  development 
the  presence  of  this  isopycnic  layer  is  accepted.  The  value  for  density 
at  90  Km  is  set  equal  to  3.44  (10-^)  gm.cm-^  for  all  months  and  latitudes 

Table  C2  depicts  how  the  January,  July  and  Spring/Fall  data  for  30°N, 

45°N  and  60°N  latitudes  are  applied  to  other  months  and  latitudes.  The 
choices  are  somewhat  arbitrary,  but  necessary  due  to  a lack  of  additional 
data. 

The  extension  of  the  density  to  400,000  feet  is  now  direct.  The 
density  at  100,000  feet  and  Tabled  are  used  to  extend  the  density 
profile  to  the  isopycnic  value  at  90  Km,  and  then  upward  to  400,000  feet. 

The  temperature  data  are  extracted  directly  from  the  U.S.  Standard 
Atmospheric  Supplements  (1966)  at  20,000  feet  increments  (Table  C3 ) • 

The  temperature  and  density  profiles  join  smoothly  with  the  Jacchia 
Model  at  120  Km.  Temperatures  and  densities  above  120  Km  can  be 
determined  from  the  Jacchia  Model  using  three  parameters:  (1)  the 
10.7-cm  solar  radio  flux  averaged  over  three  solar  rotations,  (2)  the 
daily  average  of  the  10.7-cm  solar  radio  flux,  and  (3)  the  applicable 
planetary  magnetic  index,  Ap.  All  three  of  these  solar  parameters  are 
provided  in  the  Point  Analysis. 
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A graph  of  three  actual  mid-latitude  extended  PAs  is  presented  in  Figure 
C3.  Similar  graphs  for  extended  density  profiles  are  presented  in 
Figures  C4,  C5  and  C6.  Note  that  both  the  extended  temperature  and  density 
profiles  join  the  PA  smoothly  at  100K  feet. 


TABLE  Cl 


MONTHLY  MEAN  SCALE  HEIGHTS  (KM) 

JANUARY 

30N 

45N 

60N 

30-45km 

6.810 

6.448 

6.383 

45-70km 

6.944 

7.797 

7.542 

70-8Ckm 

6.806 

7.032 

7.501 

90- 104km 

5.882 

5.979 

6.079 

104-122km 

7.358 

7.162 

7.020 

JULY 

30-4 5km 

6.894 

6.937 

6.991 

45-7Ckm 

8.058 

8.189 

8.298 

70-80km 

6.701 

6.740 

6.687 

90-100km 

5.178 

4.959 

4.862 

100-108km 

5.459 

5.420 

5.375 

108-122km 

6.800 

7.091 

7.302 

spring/fall 

30-4 5km 

6.706 

45-70km 

8.031 

70-80km 

6.733 

90-1  OCk 

5.678 

100-122km 

6.281 

ALL  MONTHS 

15N 

30-45km 

6.922 

45-7 Okm 

7.999 

70-80km  6.738 

90-108km  5.384 

108-122km  6.346 
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TABLE  C3 

MONTHLY  MEAN  TEMPERATURES  (°C) 
JANUARY 


SPRING/FALL 


pFeet) 

30N 

4^N 

60N 

120 

-29.4 

-40.1 

-46.7 

-32.3 

140 

-15.1 

-21.8 

-31.9 

-15.4 

160 

- 4.0 

- 7.5 

-16.8 

- 2.5 

180 

-10.6 

-12.3 

-13.9 

- 7.3 

200 

-24.0 

-24.3 

-22.8 

-15.5 

220 

-42.5 

-36.5 

-25.7 

-41.9 

240 

-61.0 

-49.0 

-34.5 

-62.5 

260 

-79.4 

-61.5 

-47.6 

-8O.3 

280 

-82.1 

-73.4 

-60.8 

-82.6 

300 

-75.4 

-72.6 

-68.5 

-82.0 

320 

-57.3 

-59.5 

-58.2 

-73.5 

340 

-37.0 

-41.7 

-44. 3 

-61.4 

360 

- 7.7 

-14.8 

-21.3 

-40.5 

380 

+30.0 

+26.5 

+22.7 

+22.7 

400 

+74.3 

+75.9 

+77.6 

+109.0 

TABLE  c3  (contd) 


JULY  ALL  MONTHS 


(lCTFeet) 

30N 

45N 

60N 

15N 

120 

-26.6 

-24.1 

-21.2 

-27 

140 

-12.2 

- 9.0 

- 3.2 

-13 

160 

- 1.0 

+ 2.5 

+ 3.9 

- 3 

180 

- 7.6 

- 3.5 

+ 1.1 

-10 

200 

-21.9 

-18.5 

-14.9 

-23 

220 

-44.6 

-42.1 

-43.0 

-44 

240 

-67.2 

-68.9 

-71.0 

-65 

260 

-89.8 

-95.7 

-99.1 

-86 

280 

-100.7 

-108.7 

-111.5 

-96 

300 

-98.4 

-104.6 

-107.3 

-95 

320 

-87.1 

-88.5 

-89.3 

-85 

340 

-66.7 

-62.3 

-59.4 

-69 

360 

-23.3 

-12.9 

- 5.4 

-33 

380 

+47.0 

+54.0 

+58.9 

+40 

400 

+134.0 

+130.5 

+128.5 

+138 

L 


(height  Km) 

FIGURE  C2 : Density  profile  generated  using  the  following  constant  scale 
heights  - ?,032Km  from  70  to  80  Km  - 5,979Km  from  80  to  104  Km  - 
7.l62Km  from  104  to  120  Km. 


(height  lCr1 


figure  C3 : Three  sample  extended  temperature  profiles  from  midlatitude 
Point  Analyses. 
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APPENDIX  D 


INTERPRETING 

SITE  PSEUDO-SURFACE  OBSERVATION 
(Paragraph  D.) 


D.l  The  First  Record  of  Site  Pseudo-Surface  Observation 


The  first  nonblank  record  of  paragraph  D.  contains  visibility  and  cloud 
data.  Visibility  is  interpolated  from  surrounding  surface  data  (inter- 
polation technique  is  described  in  Appendix  A) . The  cloud  data  is  taken 
from  the  nearest  3DNEPH  cloud  analysis  grid  point.  Table  D1  is  a column 
by  column  description  of  the  first  record. 

Unless  there  is  no  surface  data,  the  visibility  will  range  from  00.00  to 
15.00  nm.  When  there  is  no  surface  data  available,  99.99  is  encoded. 

The  cloud  data  is  broken  into  low,  middle,  and  high  cloud  groups.  The 
cloud  coverage,  types,  and  heights  are  presented  in  each  group.  Cloud 
coverage  is  first  and  is  given  in  eighths.  For  example,  "2/8."  Eight 
columns  are  reserved  for  cloud  types.  In  these  eight  columns,  three  two- 
character  types  are  separated  by  the  character  "/"  Table  D2  presents  all 
the  possible  cloud  types  that  may  appear.  An  example  cloud  type  code 
for  low  clouds  is  "CU/SC/**".  This  example  means  there  are  cumulus  and 
stratocumulus  clouds  present.  Seven  columns  are  used  to  give  minimum 
bases  and  maximum  tops  of  clouds  in  a layer.  Three  digits  are  used  to 
give  bases  in  hundreds  of  feet.  These  digits  are  followed  by  the  character 
"/"  and  three  more  digits  for  tops  of  clouds  in  hundreds  of  feet.  The 
characters  "***/***"  mean  missing  height  data  or  no  clouds  in  layer. 


D. 2 The  Second  Record  of  the  Site  Pseudo-Surface  Observation 

The  second  record  of  paragraph  D.  contains  wind  and  present  weather  data. 
Wind  data  in  interpolated  from  the  surrounding  surface  observations 
(interpolation  technique  is  described  in  Appendix  A.).  The  present 
weather  is  taken  from  the  closest  surface  observation.  Table  D3  is  a 
column  by  -olumn  description  of  the  second  record. 

Wind  data  is  presented  in  10  columns.  The  first  three  digits  give 
direction  in  degrees.  These  are  immediately  followed  by  three  digits 
for  speed  in  knots.  Missing  wind  data  is  encoded  as  "999999".  Calm 
winds  are  encoded  as  "000000".  These  six  columns  are  followed  by  the 
character  "G"  and  three  digits  for  gust  speed  in  knots.  Missing  gust  data 
is  encoded  as  "***".  (Gust  data  is  available  only  for  North  America.) 

Zero  to  four  12-character  present  weather  contractions  may  also  appear 
in  the  second  record.  Table  D4  presents  all  the  possible  contractions 
and  their  equivalent  in  World  Meteorological  Organization  (WMO)  code 
4677  WW.  A complete  definition  may  be  found  by  referencing  Air  Weather 
Service  Manual  105-24,  Vol.  III.  The  12-character  contractions  are 
followed  by  commas. 
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Table  Dl.  Format  for  first  Data  Record  of  Point  Analysis 


Paragraph  D. 

-cols  1-3:  three  blanks 

-cols  4-7:  the  four  characters  "VSBY" 

-col  8 : a blank 

-cols  9-13:  a four  digit  decimal  number  with  the  decimal  point  (period) 
in  col  11 

-col  14  : a blank 

-cols  15-16:  the  two  characters  "NM" 

-col  17  : a blank 

-col  18  : one  digit  for  eights  of  low  cloud  coverage. 

-cols  19-20:  the  two  characters  8 " 

-col  21  : one  blank 

-cols  22-29:  Three  two-letter  cloud  type  codes  separated  the  character 
"/" . Columns  22-29  contain  types  for  low  clouds. 

-col  30  : blank 

-cols  31-33:  three  digit  number  for  minimum  low  cloud  bases  in  hundreds 
of  feet  above  ground  level. 

-col  34  : the  character  "/" 

-cols  35-37:  three  digit  number  for  the  maximum  low  cloud  tops. 

-col  38  : a comma 

-cols  39-58:  same  as  columns  18-37  except  that  the  data  refers  to 
middle  clouds 

-col  59  : a comma 

-cols  60-79:  same  as  column  18-37,  except  that  data  refers  to  high 
clouds 

-col  80  : blank 
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Table  D2.  Possible  cloud  types 


TyP-e- 


SC 

ST 

CU 

CB 

AC 

AS 

NS 

Cl 

CC 

CS 

** 


UN 


Meaning 

Stratocumulus 

Stratus 

Cumulus 

Cumulonimbus 

Altocumulus 

Altostratus 

Nimbostratus 

Cirrus 

Cirrocumulus 
Cirrostratus 
Missing  data, 

or  no  other  cloud  type  in  layer 
or  no  clouds  at  all  in  layer 

Unknown  (the  3DNEPH  cloud  analysis 
program  was  not  able  to  determine 
the  cloud  type) 
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Table  D3.  Format  for  the  Second  Data  Record  of 
Paragraph  D.  in  the  Point  Analysis  Product 

cols  1-3:  three  blanks 

cols  4-7:  the  four  characters  "WIND" 

•col  8 : a blank 

cols  9-11:  three  digits  giving  wind  direction  in  degrees, 
cols  12-14:  three  digits  giving  wind  speed  in  knots 
col  15  : the  character  "G" 

cols  16-18:  three  digits  giving  gust  speed  in  knots 
cols  19-21:  three  blanks 

cols  22-73:  one  to  four  12-character  present  weather  contractions 
each  one  followed  by  a comma 

cols  74-80:  seven  blanks 


..  b£ST  quality  mCXlCABLB 

THIS  PAGE  Is  BE  jo  ddc  

7R0M  COP*  FUBIilSHEO  ^ 

Table  D4.  List  of  All  Possible  Twelve  Character 
Present  Weather  Contractions  and  their 
Equivalent  in  WMO  Code  4677  WW. 

(Reference  AWSM  105-24,  Vol.  Ill) 


COM  TRACT  ION 

WHO 

*677 

— m 

on 

Cl o ncRSC-  - 

0 1 

- cld  same 

U 2 

“ C L 0 BLOG 

U3 

5M0KL 

U 4 

HAZL 

US 

- OORT 

0 6 

~ BL  wG  OUST 

U 7 

- dost  dev i i s- 

U8 

- DOST  STORM  - 

U 9 

- MIST 

i n 

“PC  M Y GRN  F'Og- 

1 1 

“ GROUND  FOG  - 

l 2 

- I TG  5tFN 

i 3 

vIRGa 

1 4 

-DISTANT  PRCP- 

1 S 

- CLOSE  PRCP  - 

1 6 

-THMoR  HEARD  - 

1 7 

- SwUALLS 

1 8 

- FUNNEL  cld  - 

1 R 

-DR7l/SNw  GRN- 

20 

“ R A I N 

2 i 

- SNO* 

22 

“RN/SNW/ICEPL- 

23 

- FPZG  PRcP  - 

24 

- rain  shhrs  - 

2 S 

-rn/snw  shwr  - 

2 6 

-RN/HAIL  n X 0 - 

27 

“ FDg/ICE  FOfe- 

2R 

-ThUnDERSTORM- 

29 

-DCRSG  MOOST  - 

JO 

- moo  oust 

31 

-incrsg  mdust- 

3? 

-DCRSG  SOOST  - 

33 

- s v r dost 

34 

-INCRSG  5DUST- 

35 

-md  uRFTo  SNk- 

36 

-HV  oRFTg  SNw- 

37 

-M0D  BUWg  snw- 

3 P 

•HVY  BLW(j  SNw- 

3 9 

■ DS  TNT  FOG  - 

40 

- PTtHY  FOG  - 

4 1 

- Stfr/pc  FOG  - 

4? 

■N0SkY/0c  FOo- 

4 3 

■ sky  fog 

4 4 

■ NO^KY  fog 

4 S 

•sky/in  fog 

46 

■NOSkY/In  F0(,. 

47 

• SKY/RM  FOG  - 

4 R 

-NOSnY/RFi  FDg- 

49 

CONTRACT  I ON 

WHO 

46/7 

-L'jT  INT  DR7L- 

SO 

-LOT  CON  DRZL- 

5 I 

-Mur  TNT  0R7L- 

52 

-Mor  CON  0R7L- 

6 3 

-HVY  TNT  ORZL- 

54 

•HVY  CON  Ph  7L- 

55 

-LGT  F R 7 PR7L- 

66 

-MOD  F R 7 OR  ZL- 

67 

■LT  DRZL+RaIN- 

5 8 

■MD  DRZL*RAiN- 

69 

■lgt  InT  rajn- 

«n 

-lgt  con  rain- 

6 1 

•MOD  JNT  RA1N- 

62 

-iion  con  rain- 

63 

*H  V V InT  RaIN- 

64 

■HVY  foN  RaIN- 

65 

■lgt  F R 7 RaIN- 

66 

•MOD  F R 7 RaIN- 

67 

■L  T KaIN+SNCW- 

6 R 

■nd  R A I N ♦ S N 0 W - 

69 

■L(,T  |NT  SNOW- 

70 

■LGT  CON  snow- 

7 1 

■FiOO  InT  SNOW- 

72 

■MOD  CON  SNOW- 

/ 3 

•HVY  |fjT  SNOW- 

7 4 

■HVY  CON  SNOW- 

75 

ICE  CRTSTALS- 

76 

■SNOr  f,RAINS  - 

77 

■ISLTp  CRYSTL- 

7 ft 

• I CE  PELLETS  - 

79 

■i.T  rain  shwr- 

bO 

•mdh  r a i n shw- 

8 l 

■S  V R RaIn  shw- 

82 

LT  RAIN+SNOW- 

83 

MD  KA|N*SNOW- 

84 

LT  SNOW  ShWR- 

tt6 

NDH  SNW  SHWR- 

86 

LT  SNOW  PLT  - 

87 

MDH  SNOr  Pt  T- 

88 

LT  HAIL  SHWR- 

89 

P'DH  HaIL  shw- 

VO 

LT  rain  nbT  - 

9 1 

NDH  RaIn  GpT- 

V 2 

LT  SNOW  0 b T - 

9 3 

MDH  snow  OpT- 

V 4 

T H D S T M ♦ R A J N - 

96 

1 nDSTM*KN«.HL- 

46 

H V Y ThDSTN  - 

V 7 

THOSTM^DUST  - 

98 

HV  1HdSTM*hL- 

99 
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APPENDIX  E 

INTERPRETING  THE 
SURFACE  OBSERVATIONS  TABLE 
(Paragraph  G.) 


67 


r 


E.l  Inclusion  of  Surface  Observation  Table 


The  PA  user  may  request  that  data  from  surface  observation  stations 
surrounding  the  input  site  be  included  in  Paragraph  G.  The  stations 
chosen  are  as  close  as  possible  to  the  user's  input  site  and  time. 

E.2  Table  Entries 


All  possible  types  of  data  are  listed  in  TableEl.  Each  line  of  the 
observation  table  in  paragraph  G.  is  labeled  as  shown  in  Table  El. 
Entries  for  which  no  data  is  available  are  deleted. 

The  entries  for  "BLOCK/STATION  NUMBER"  through  "DEWPOINT  (DEG  C) " 
are  self-explanatory.  The  user  may  request  inclusion  of  the  surface 
observation  table  with  the  entries  for  "BLOCK/STATION  NUMBER"  through 
"BEARING  (DEG)"  deleted. 


E.2.1  Present  Weather 


Up  to  four  present  weather  contractions  may  appear  for  each  station. 
The  contractions  are  exactly  the  same  as  those  shown  in  Appendix  D, 
Table  D4 


E.2. 2 Past  Weather 


One  twelve-character  past  weather  contraction  may  appear  for  each 
station.  Table E2  lists  all  the  possible  contractions  and  their 
equivalent  in  WMO  code  4500  (reference  AWSM  105-24,  Vol.  Ill) . 


E.2. 3 Clouds 


There  are  four  possible  kinds  of  entries  under  the  data  type  "CLOUDS" 

The  first  type  is  the  principal  cloud  type  the  observer  reported  for 
low,  middle  and  high  clouds.  The  distinguishing  trait  of  this  type 
is  that  two  commas  appear  in  it.  Two  characters  and  a digit  are  encoded 
for  each  layer.  The  single  digit  is  from  WMO  code  513  for  low  clouds, 

WMO  code  515  for  middle  clouds  and  WMO  code  509  for  high  clouds 
(reference  AWSM  105-24,  Vol.  III).  The  characters  "MIS"  indicate 
missing  data. 

The  second  type  of  "CLOUDS"  entry  gives  cloud  coverage,  cloud  type, 
height  of  cloud  bases  for  important  cloud  layers  reported  by  the  observer. 
The  distinguishing  trait  of  this  type  is  the  presence  of  one  comma. 

There  may  be  zero  to  six  entries  of  this  type  for  each  station.  The 
cloud  coverage  may  vary  from  zero  to  eight  eights.  If  the  sky  is 
obscured  "9/8"  will  appear.  The  cloud  type  code  is  identical  to  that 
used  above.  Three  digits  give  cloud  bases  in  hundreds  of  feet  above 
ground  level. 


68 


J 


The  third  type  of  "CLOUDS"  entry  gives  the  total  fraction  of  the  sky 
covered  by  clouds.  Its  distinguishing  trait  is  the  characters  "TOTL 
CLD".  The  cloud  coverage  is  given  in  eighths.  If  the  cloud  ceiling  is 
indefinite  "9/8"  will  appear.  If  the  observer  reported  a sky  obscura- 
tion, the  characters  "SKY  OBSCURED"  will  replace  this  entry. 

The  last  type  of  "CLOUDS"  entry  gives  the  lowest  reported  ceiling. 

Its  distinguishing  trait  is  that  it  begins  with  the  characters  "CIG". 

The  height  of  the  lowest  ceiling  is  given  in  hundreds.  The  three 
digits  for  height  may  be  preceded  by  a single  letter  which  either 
indicates  how  the  ceiling  height  was  determined  (M-measured,  E-estimated, 
A-aircraft,  or  B-ballon) , or  characterizes  the  ceiling  (W-indef inite, 
X-obscuration) . If  the  ceiling  is  "variable"  or  "thin",  then  the  last 
four  characters  of  the  entry  will  be  "VRBL"  or  "THIN",  respectively. 
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TABLE  El 


LIST  OF  ALL  POSSIBLE  DATA  TYPES 
IN  THE  SURFACE  OBSERVATIONS  TABLE  OF  PARAGRAPH  G. 


BLOCK/STATION  NUMBER 
CALL  LETTERS 
ZULU  TIME 
BEARING  (DEG) 

DISTANCE  (NM) 

WIND  DIRECTION  (DEG) 
WIND  SPEED  (KTS) 

WIND  GUSTS  (KTS) 

SEA  LEVEL  PRESSURE  (MB) 
3-HRLY  PRES  CHANGE  (MB) 
ALTIMETER  SETTING  (IN) 
VISIBILITY  (NM) 

6-HR  PRECIPITATION  (IN) 
TEMPERATURE  (DEG  C) 
DEWPOINT  (DEG  C) 

PRESENT  WEATHER 
PAST  WEATHER 
CLOUDS 
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TABLE  E2 


LIST  OF  ALL  POSSIBLE  TWELVE  CHARACTER  PAST  WEATHER 
CONTRACTIONS  AND  THEIR  EQUIVALENT  IN  WMO  CODE  4500. 
(REF  AWSM  105-24,  VOL  III) 


CONTRACTION  tr  M 0 45  0 0 


-LT  1/?  CLOY  - U 

- Y R T i 

-gT  t/?  f LD  Y - ? 

- DU5T5TORM  - 3 

- FOG  - 4 

- 0RI7ZLL  - b 

HAJN  - 6 

- bNOW  - 7 

SH0WEP5  - H 

-IHUNUFPSTOHM-  9 
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APPENDIX  F. 

USER  CONSIDERATIONS 


F . 1 Protecting  Input  Site  and  Time 

The  basic  Point  Analysis  message  format  clearly  reveals  the  user’s 
input  site  and  time.  If  the  user  wants  the  site  and  time  protected,  then 
the  following  changes  are  made: 

Para.  A.  - Site  latitude  and  longitude  are  replaced  by  a six 
character  identification  code  provided  by  the  user  at  the  time  of  his 
request. 

Para.  B.  - Paragraph  B.  is  left  blank. 

Para.  G.  - If  a rawinsonde  was  used  to  create  the  vertical  profile 
the  rawinsonde  station  number  will  be  deleted. 

Para.  G.  - If  the  user  has  requested  inclusion  of  the  Surface 
Observations  Table,  the  entries  for  BLOCK/STATION  NUMBER  through  "BEARING  (DEG)" 
are  deleted. 

F . 2 Analysis  Only  Versus  Rawinsonde 

The  users  of  Point  Analysis  products  desire  a vertical  profile  that 
is  as  representative  as  possible.  For  some  applications,  and  for  some 
areas  of  the  world  where  data  is  sparse,  it  is  advisable  to  force  the 
PA  program  to  ignore  rawinsonde  data  and  build  the  vertical  profile  based 
on  data  at  the  closest  GWC  analysis  fields  grid  point. 

However,  the  GWC  Analysis  Fields  used  by  the  PA  program  are  the 
worldwide  whole  mesh  fields.  The  grid  spacing  for  these  fields  is 
approximately  200  nautical  miles  at  60°  latitude.  Therefore,  the  closest 
grid  point  can  be  a maximum  of  150  nm  away,  and  average  about  70  nm 
away.  Furthermore,  only  the  standard  pressure  levels  have  data  and  never 
more  than  four  levels  have  moisture  data  (see  Figure  6 in  manual) . 

Most  PA  runs  find  a rawinsonde  with  at  least  four  levels  of  moisture 
data  within  150  nm  of  the  input  site.  Therefore,  allowing  the  PA  program 
to  search  for  Rawinsonde  before  falling  back  to  Analysis  Fields  data  is 
recommended.  This  is  especially  so  for  applications  which  require  the 
vertical  stratification  of  density  and  moisture  available  in  the  Rawinsonde 
significant  levels  data. 

F. 3 Summary  of  Forecaster  Capabilities 

Because  some  human  judgement  of  what  is  or  is  not  representative  of 
a given  site  and  time  is  useful,  the  PA  production  process  was  designed 
to  allow  GWC  forecaster  to  intervene  in  the  process. 
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F . 3 . 1 Selection  of  Data  for  the 

The  forecaster  can  force  the  PA  program  to  ignore  1 to  5 rawinsondes 
which  he  feels  are  not  representative.  He  can  also  look  at  a listing  of 
available  rawinsondes  and  force  the  PA  program  to  use  a particular  rawinsonde 
which  the  forecaster  feels  is  representative.  Frequently,  the  program 
determines  there  is  no  rawinsonde  available  and  falls  back  to  GWC  Analysis 
Fields  data  on  its  own.  However,  the  forecaster  can  also  decide  to  ignore 
rawinsonde  data  and  force  the  program  to  do  so. 

F . 3 . 2 Change  the  Pseudo-Surf ace  Observation 

The  GWC  forecaster  may  review  the  contents  of  paragraph  D.  and  determine 
some  part  of,  or  all  of,  the  Pseudo- surface  Observation  not  representative. 

He  may  change  any  part  of  this  paragraph  to  whatever  he  feels  is  representative. 

F. 3. 3 Paragragh  G Remarks 

The  GWC  forecaster  is  free  to  include  whatever  plain  language  remarks  he 
desires.  This  will  include  any  special  information  requested  by  the  user, 
information  which  the  forecaster  feels  the  user  ought  to  be  aware  of,  or 
general  remarks  about  the  representativeness  of  the  vertical  profile. 
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